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We consider interacting electrons in a two-dimensional quantum Coulomb glass and investigate by 
means of the Hartree-Fock approximation the combined effects of the electron-electron interaction 
and the transverse magnetic field on fluctuations of the inverse compressibility. Preceding systematic 
study of the system in the absence of the magnetic field identifies the source of the fluctuations, 
interplay of disorder and interaction, and effects of hopping. Revealed in sufficiently clean samples 
with strong interactions is an unusual right-biased distribution of the inverse compressibility, which 
is neither of the Gaussian nor of the Wigner-Dyson type. While in most cases weak magnetic 
fields tend to suppress fluctuations, in relatively clean samples with weak interactions fluctuations 
are found to grow with the magnetic field. This is attributed to the localization properties of the 
electron states, which may be measured by the participation ratio and the inverse participation 
number. It is also observed that at the frustration where the Fermi level is degenerate, localization 
or modulation of electrons is enhanced, raising fluctuations. Strong frustration in general suppresses 
effects of the interaction on the inverse compressibility and on the configuration of electrons. 

PACS numbers: 73.23.Hk, 71.55.Jv, 71.30.+h 



I. INTRODUCTION 

Since the experimental observation of the crucial ef- 
fects of the Coulomb interaction,El mesoscopic fluctua- 
tions of the inverse compressibility in two-dimensional 
quantum dots have been i nv estis^(^|iA^Htti]*G|ly, both 
theoretically., and experimentally.l!lal3'QBl3t3'ErBli3 Recent 
experimentaj'BEl on the distribution of the inverse com- 
pressibility, which can be measured by conductance peak 
spacing, in the Coulomb-blockade regime have shown 
that the standard random matrix theory fails to ex- 
plain the observed fluctuations; this implies that charg- 
ing energy fluctuations may play a dominant role in the 
inverse-compressibility fluctuations. In particular, the 
distribution was observed to take a|-jGaussian-like sym- 
metric form with non-Gaussian tailsMacI Subsequent an- 
alytical and numerical studies on interacting electrons 
in disordered systems revealed the possibility of larger 
inverse-comOT|es.silwlity fluctuations induced by Coulomb 
interactions BqBQID'EI It was also proposed that the shapg 
deformation of the dot due to the-gate- voltage sweepings 
or the irregular shape of the dottL3 may contribute to the 
deviation of the experimental results from the Wigner- 
Dyson (WD) statistics predicted by the random matrix 
theory. However, the recent experiment with a stacked 
gate structure, which allows one to vary the electron den- 
sity without significant deformation of the dot, was sup- 
portive of the claim that the observed shape of fluctua- 
tions originates from electron-electron interactionsfl An- 
other recent experiment reported the crossover behavior 
of the compressibility axound the metal-insulator transi- 
tion in two dimensions JIj 

The influeiice of the .magnetic field, either 
perpendicular™!! or parallel)! 2 ! on the distribution 



of the inverse compressibility was also investigated 
experimentally. In the presence of a perpendicular 
magnetic field, the shape of the distribution is found 
qualitatively similar to that in zero magnetic field, 
only with narrower width. u There was an attempt to 
explain the qualitative behavior of psak positions in 
the system under the magnetic fieldjlij the considered 
classical model is based on the electrostatics of several 
electron islands with quantum orbitals, the energies 
of which depend on the external magnetic field. The 
possibility of nearly vanishing peak spacing was thus 
suggested, which is in qualitative .-agreement with the 
bunching of the addition spectra.E-3 However, it does 
not take into account the quantum interference effects 
between electron wave functions, which can be crucial in 
mesoscopic systems. 

In this paper, we study numerically fluctuations of the 
inverse, p^mpressibility in the two-dimensional Coulomb 
glass. EjEEIEZI To treat the interaction between electrons, 
we employ the Hartree-Fock (HF) approximation, which 
makes it possible to investigate larger samples than those 
feasible in exact diagonalization. For a systematic study, 
we first examine the classical system, revealing the de- 
tailed dependence of the inverse-compressibility fluctu- 
ations upon disorder and interaction. Then the ef- 
fects of electron hopping on the fluctuations are inves- 
tigated. Observed is an unusual right-biased distribution 
for very weak disorder and strong interactions as well as a 
symmetric Gaussian-like distribution with non-Gaussian 
tails. It is also found that relative fluctuations decrease 
eventually for sufficiently strong interactions, which ap- 
pears to be inconsistent with the quadratic interaction 
dependence of the fluctuations suggested in a previous 
numerical studyfl We further examine the effects of mag- 
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netic fields, which introduce frustration to the system, 
on the distributions and fluctuations of the inverse com- 
pressibility, and observe that weak magnetic fields tend 
to suppress fluctuations. Nevertheless the fluctuations 
may be enhanced in some range of the frustration for 
weak disorder and interaction. Such opposite tendency 
of the inverse-compressibility fluctuations in response to 
the applied magnetic fields is explained in terms of the 
localization properties of the system. Finally, it is also 
found that strong frustration suppresses significantly the 
effects of the interaction. 

This paper is organized as follows: In Sec. II we con- 
sider the classical Coulomb glass and investigate in de- 
tail the effects of disorder and interaction on the inverse- 
compressibility fluctuations. The effects of hopping on 
the fluctuations in the quantum Coulomb glass are ex- 
amined in Sec. [II. Section IV is devoted to the effects 



of magnetic fields. We examine the system at both weak 
and strong frustration, and describe the electron distribu- 
tions and localization properties by means of the partic- 
ipation ratio and the inverse participation number. The 
main results are summarized in Sec. 



II. CLASSICAL COULOMB GLASS 

We address,the Hamiltonian for the classical Coulomb 
glass modelfj 



\Ytfk-K)U ij {b i -K), (1) 



where i is the site index on an L x L square lattice, hi 
is the number operator of the electron at site i, and Wi 
is the random on-site potential uniformly distributed in 
the range [-W, W]. We assume the Coulomb interaction 
potential between electrons: 
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with e being the dielectric constant. In Eq. ([j]) K rep- 
resents the uniform neutralizing background charge per 
site in units of +e, which is demanded by the charge neu- 
trality condition over the whole system. Throughout this 
paper we concentrate on the half-filling case K = 1/2 in 
the free boundary conditions, and write the distance in 
units of the lattice constant a. 

The ground state of the classical Hamiltonian in 
Eq. (Q) is obtained by means of simulated annealing pro- 
cedure. The inverse compressibility A2 for a given disor- 
der configuration is then calculated via the relation 
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where M = L 2 /2 and E n denotes the n-electron ground- 
state energy. We have performed the calculation for size 
L = 4, 6, and 8, and averaged the data over 2000 different 
disorder configurations. 
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FIG. 1: (a) Averages and (b) fluctuations of the inverse com- 
pressibility in units of V as functions of W/V in the classical 
system of size L — 4(D), 6(0)i an d 8(A). Each straight line 
represents the behavior of the corresponding noninteracting 
system (V = 0). The insets display the averages and the 
fluctuations in units of W as functions of V/W. 



Figure [l](a) shows the average inverse compressibility 
as a function of the relative disorder strength W/V. For 
weak disorder, (A2) /V depends rather weakly on both 
the disorder strength W and the size L since (A 2 ) is pro- 
portional to V in this regime as shown in the inset. As 
the disorder is increased, (A2) /V is suppressed up to 
W/V ~ 1, then increases for W/V > 1, and approaches 
asymptotically the straight line describing the behavior 
in the corresponding noninteracting system. One can 
also observe that (A2) decreases with the system size, 
which is prominent for strong disorder. It is remarkable 
that the disorder-scaled average (A2) /W, which is plot- 
ted in the inset of Fig. |l|(a), is monotonically increasing 
with V/W , thus indicating that larger interactions lead 
to larger values of the inverse compressibility for given 
disorder strength. 

In contrast to the average (A 2 ), the fluctuations SA 2 
of the inverse compressibility, shown in Fig. |l|(b), scale 
with W independently of V for weak disorder as well as 
for strong disorder. In the crossover regime, on the other 
hand, the fluctuations display non-monotonic behavior, 
with a broad peak at W/V w 0.2. This becomes conspic- 
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FIG. 2: Relative fluctuations versus W/V in the classical sys- 
tem. The symbols are the same as those in Fig. |l|. 



uous as the system size grows: Whereas for W/V > 0.2, 
SA 2 decreases with the size similarly to (A 2 ), slight in- 
crease of SA 2 can be observed for W/V < 0.2. It is of 
interest to note that SA 2 approaches the noninteracting 
limit rather fast as W/V is increased, demonstrating that 
for strong disorder (W/V ^> 1) the interaction merely 
shifts the average inverse compressibility toward slightly 
larger values. 

In Fig. U we plot the relative fluctuations SA 2 / (A 2 ) 
of the inverse compressibility versus the relative disorder 
strength W/V. The relative fluctuations in general grow 
with W/V, saturating to a finite value in the noninteract- 
ing limit. At the point W/V = 1 the relative fluctuations 
are apparently independent of the system size_which has 
been reported in previous numerical studiesBQ However, 
it should be noted that such universal behavior shows up 
only at W/V = 1: While the relative fluctuations grow 
with the system size for W/V < I , they are reduced by 
the increase of the system size for W/V > f. Again 
a broad peak, which becomes prominent for larger sys- 
tems, is observed for W/V ss 0.2; this can be attributed 
to the presence of the peak in <5A2 mentioned above. 

In the limit of strong disorder (W/V ^> f), the elec- 
trons are distributed spatially to minimize the total ran- 
dom potential energy, and firmly pinned, rarely rearrange 
their whole distribution in response to the addition or 
removal of an electron. Consequently, as in the nonin- 
teracting system, both (A 2 ) and 5A 2 are expected to 
increase with the disorder strength W, while the interac- 
tion simply renormalizes the values slightly without qual- 
itative change. Note that such tendency begins at rela- 
tively small disorder W/V > 0(1), as shown in Fig. [|. 

In the opposite limit of strong interaction (W/V <C 1), 
on the other hand, the half-filled electron state \M) is 
expected to form a Wigner crystal (WC). In case that 
two WC states are nearly degenerate even in the presence 
of the disorder, the addition or removal of an electron can 
change the ground state into that similar to the other WC 
state; however, the disorder configuration, which allows 



such phenomena, is very rare and accordingly, the WC 
state in general does not experience significant change 
by the insertion or extraction of an electron. The inverse 
compressibility can then be calculated as follows: 

A 2 = (Em+i — Em) + (Em -i — Em) 

= w l+ + y^X n i ~ K )Uu + - w t _ - y^(rtj - K)U vi _ 



= w i+ - Wj_ + }Xnj - K)(U ii+ - Ug_), (4) 
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where rij is 0/1 for site i empty/occupied in the state 
\M), and i± represents the site to/from which an elec- 
tron is inserted/extracted. In this limit, the interaction 
part plays a major role in choosing i±, which is very likely 
to be one of the corners. Accordingly, the strength w, ± 
of the impurity at site i± is random and uncorrelated, 
leading to the disorder-independent average given by the 
third term on the right-hand side of Eq. (|4|) ; this quantity 
decreases slightly with L. In contrast, the fluctuations 
arise mainly from the first two terms which have their 
origin in the random potential and thus increase with 
the disorder strength W. The numerical data presented 
above are indeed in good agreement with these conclu- 
sions based on analytic argument. Figure [l] demonstrates 
that this regime persists until W/V w 0.02. 

In both extreme limits considered above, one can see 
that the contributions of the random potential to the 
fluctuations of the inverse compressibility are dominant 
over those of the interaction. In the intermediate re- 
gion, however, the behavior looks rather complicated and 
the interaction is expected to play an important role. 
For detailed investigation, we divide the intermediate re- 
gion into three subrcgions: (I) 0.02 < W/V < 0.2, (II) 
0.2 < W/V < 0.5, and (III) 0.5 < W/V < 0(1). In Re- 
gion I the half- filled state \M) strongly depends on the 
disorder configuration: Some states are still close to the 
WC state while others are rather disordered. In the lat- 
ter states, the portion of which increases with W/V, the 
distribution of electrons is inhomogeneous, reducing the 
Coulomb energy cost of the addition or removal of one 
electron. On the other hand, the former states, which 
resemble the WC state, rather easily rearrange the dis- 
tribution of electrons by the addition or removal of an 
electron, unlike in the limit of strong interactions. Both 
behaviors generally result in the reduction of A 2 although 
the amount of the reduction depends upon the disorder 
configuration; this explains the change of the A 2 /V dis- 
tribution from the 5-function-like one [see Fig. [|(a)] to 
the asymmetric broad distribution with a long left tail 
[see Fig. ||(b)]. It is observed that the minimum of A 2 /V 
as well as the peak position of the distribution decreases 
with W/V. Indeed in a finite system the minimum of 
A 2 /V hits its lower bound given by the smallest pos- 
sible Coulomb interaction between successively inserted 
electrons ,Q 
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FIG. 3: Distribution of A2/V in the classical system of size L = 8 for disorder strength W/V = (a) 0.01, (b) 0.1, (c) 0.2, (d) 
0.3, (e) 1, (f) 5. 



at W/V ~ 0.2. In addition, as W/V is raised, the 
emergence of rather disordered states in turn leads to 
the increase of the fluctuations of A2/V, which is re- 
flected as well in the correlations between Em+i — Em 
and Em-i—Em shown in Fig. ^. In this regime the corre- 
lations are positive and increase with W/V until W/V « 
0.2, keeping parallel with the broader distribution of 
A 2 /V. [Recall that A 2 = {E M+l -E M ) + [E M -i-E M ).\ 
Such positive correlations may be understood via the ar- 
gument that the presence of a dilute region in state |M), 
which lowers Em+i — Em, accompanies a dense region, 
allowing a low cost of the electron extraction. It should 
be noted that large fluctuations in this regime arise from 
the interplay between the Coulomb interaction and the 
random potential. 

The main feature of Region II is that essentially all the 
states are disordered. In this regime the disorder contri- 
butions become important in determining the electron 
configuration of |A/), alleviating the need for rearrange- 
ment due to the insertion or extraction of an electron. 
The fluctuations S[Em±i — Em], which increase or de- 
crease very slightly depending on the system size, are 
not enough to explain the noticeable decrease of 5A2/V 
shown in Fig. 0(b), suggesting that the correlations be- 
tween them play a crucial role. Figure ^ shows that the 
correlations decrease to negative values. This behavior, 
providing an evidence for the mounting role of the ran- 
dom potential, can be understood as follows: In the non- 
interacting case Em+i—Em and Em-\—Em are simply 
given by wm+i and —wm, respectively, where w n is the 
nth smallest strength of the on-site potential in a dis- 
ordered sample. The closeness of wm+i and wm gives 
positive correlations between them, thus leading to neg- 
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FIG. 4: Correlations between Em+i— Em and Em-\—Em in 
units of V 2 in the classical system of size L — 4(D), 6(0)i 
and 8(A). 



ative correlations between Em±i — Em- This argument 
is expected to remain valid in the presence of sufficiently 
weak interactions. It is remarkable that such correlations 
indeed begin to show up in the regime of rather strong 
interactions. Furthermore, as can be seen from Fig. |3](c) 
and (d), the distribution of A 2 /V is bounded below by 
the minimum of A 2 /V given by Eq. (|5|), while the max- 
imum of A 2 /V decreases with W/V, accounting for the 
corresponding decrease of the fluctuations SA 2 /V. 

Finally, one can observe that Region III is charac- 
terized by weak dependence of (A 2 ) /V, SA 2 /V, and 
SA 2 / (A 2 ) upon the disorder strength. Although fluc- 
tuations of Em+i—Em and of Em-\—Em increase with 
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W/V , they are canceled out due to their negative cor- 
relations as disclosed in Fig. [|. It should also be noted 
that the inverse-compressibility distribution follows the 
WD form throughout this regime [see Fig. ||(e)] while the 
universal relative fluctuations reported in Refs. |^ and |^ 
emerge only for W/V = 1. 

Before closing this section, we add a remark on the 
shape of the distribution of the inverse compressibil- 
ity, which is displayed in Fig. || for several values of 
W/V. The noninteracting system follows the Pauli dis- 
tribution, which changes to the WD distribution for 
0.5 < W/V < 0(1) as pointed out above. Further de- 
crease to 0.05 < W/V < 0.2 generates quite a peculiar 
distribution: the asymmetric one with a peak biased to 
the right and a long tail in the left, which resembles nei- 
ther the WD nor the Gaussian distribution. As W/V be- 
comes even smaller, this distribution eventually reduces 
to the (5-function distribution, manifesting the WC state. 

III. QUANTUM COULOMB GLASS 

In this section we allow hopping of electrons between 
sites and investigate the corresponding quantum effects. 
The Hamiltonian for such a quantum Coulomb glass 
readsBEl 

H = t + CjCi) + ^ W i™i 

(id) '' 

+ 1 -Y J (n i -K)U ij {n j -K), (6) 

where c\/ci is the creation/annihilation operator at site i 
and the number operator is given by hi = cjcj. The hop- 
ping of an electron is allowed only between nearest neigh- 
boring sites as indicated by the first term. We employ 
the HF approximation to treat the Coulomb interaction 
in Eq. (||), and obtain 

Hhf = ty~)(c]cj + c]cj) + ^ w iiii 

(id) 1 

apart from constant terms. The third and the fourth 
terms represent the direct and the exchange energies of 
the electron-electron interactions, respectively, and the 
angular brackets denote the expectation values with re- 
spect to the HF ground state, which should be deter- 
mined in a self-consistent manner. During the iteration 
we use the method of potential mixing, where the aver- 
ages over a few previous iterating steps are inserted into 
the expectation values in the Hamiltonian; this not only 
ensures convergence of the iteration but also helps to re- 
duce iteration steps necessary for convergence. 

We have computed the inverse compressibility up to 
2000 disorder configurations for several values of the hop- 



ping strength t/V and the disorder strength W/V. Fig- 
ure |5| displays the distribution of the inverse compressibil- 
ity for various values of t/V and W/V, together with the 
distribution in the corresponding classical system. For 
small W/V shown in Figs. ||(a), (b), and (c), as electron 
hopping comes into play, the peak position of the dis- 
tribution moves toward smaller values of A2, with the 
probability for large A2 reduced. One can see that for 
W/V = 0.1 the presence of even very weak hopping pro- 
duces a noticeable change in the distribution. Indeed 
the clean system (W — 0) is observed to show rapid 
decrease of the inverse compressibility with the increase 
of the hopping strength, suggesting that the decrease 
of the probability for large A2 originates from the in- 
terplay between the Coulomb interaction and hopping. 
Figure ||(d), (e), and (f) corresponding to W/V > 0.5, 
in contrast, demonstrates that the peak shifts slightly 
to the right and appears to saturate eventually as the 
hopping strength is increased. In addition, the probabil- 
ity for small A2 is suppressed and the minimum of A2 
rises, which appears to begin at somewhat weaker disor- 
der W/V ~ 0.2. Recalling that such suppression in the 
probability for small A2 is also observed in the noninter- 
acting system, we can infer that random disorder plays a 
crucial role here. The effects of electron hopping are also 
manifested by the statistics of the inverse compressibil- 
ity, which is presented in Fig. ra. For W/V < 0.5, both 
(A2) and 5A2 decrease with t/V due to the reduction 
of the probability for large A2; for W/V > 0.5, on the 
other hand, the suppression of small A2 induces slight 
increase of (A 2 ) and decrease of SA 2 , which results in 
overall suppression of 6A2/ (A2) in the whole region. 

Figure ^ illustrates the interaction effects on the 
inverse-compressibility distribution for given hopping 
and disorder strength. As the interaction strength is in- 
creased, the peak in the distribution of the normalized 
inverse compressibility (A2 — (A2))/<5A2 moves to the 
right both in the clean case (W/t = 1) and in the dirty 
case (W/t = 5). Here it is of interest to notice the growth 
of a tail in the left side. The increase of the interaction 
modifies the WD-like distribution to a symmetric one and 
even drives the distribution into a right-biased shape at 
sufficiently strong interactions [see Fig. 0(c)]. The distri- 
bution returns to a symmetric narrow one upon further 
increase of the interaction. Strong disorder in dirty sam- 
ples, on the other hand, tends to suppress such influence 
of the interaction: Growth of the left tail is retarded and 
the distribution remains in the left-biased shape even at 
strong interactions. 

Several experimentsSSi reported that the distribu- 
tion of the normalized inverse compressibility follows 
a roughly symmetric Gaussian form with non-Gaussian 
tails, which is very different from the WD distribution 
predicted by the random matrix theory. The interaction 
strength in those experiments is characterized by the di- 
mensionless parameter r s ss 1 to 2, which, via the relation 
r s = V/(2ty/vir) with the filling factor v = 1/2, corre- 
sponds to V/t a 2 to 5 in our model. In this interaction 
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FIG. 5: Distributions of A 2 /V in the quantum system of L = 8 for W/V = (a) 0.1, (b) 0.2, (c) 0.3, (d) 0.5, (e) 1, and (f) 
2. For comparison, the distribution in the classical system (t/V = 0) is plotted by a thick solid line in each figure. The data 
marked by squares, circles, triangles, inverted triangles, and diamonds correspond to the hopping strength t/V = 0.1,0.5, 1,2, 
and 10, respectively. 



range and also at weaker interactions we have identified 
such a Gaussian-like symmetric distribution as long as 
disorder is sufficiently weak (W/V < 0.3). In addition 
to this symmetric distribution, we have also observed 
an interesting asymmetric right-biased distribution with 
long left tails for 0.05 < W/V < 0.2 and t/V < 0.2 
[see Fig. 7(c)]. Indeed, such an asymmetric distribution, 
which also appears in Fig. 3(b) iox. the classical system, 
has been observed experimentallyo in the samples with 
strong interactions r s « 2.1 and rather weak disorder. 
[See Fig. 3 in Ref. Furthermore, we have obtained in 
this regime relative fluctuations smaller than the vak*es, 
0.1^0.2 obtained in previous numerical calculationspoQ 
this allows one to explain the weak relative fluctuation 
(5A 2 / (A2) ~ 0.06 observed in experiment^ without con- 
sidering thermal broadening. All these observations thus 
support the conclusion that the peak spacing fluctuations 
in the Coulomb blockade conductance originate mainly 
from charging energy fluctuations. 

We also investigate the effects of the interaction on the 
statistics of the inverse compressibility by computing the 
relative fluctuations, which are shown as a function of V/t 
for several values of W/t in Fig. ||. It is observed that 
particularly for intermediate disorder strengths a broad 
peak develops at a certain value of V/t, which turns out 
to increase with W/t. The decrease of <5A 2 / (A 2 ) at suf- 
ficiently strong interactions has its origin in the fast in- 
crease of (A2) as well as the decrease of 6A2 with V/t in 
this regime, where the half-filling state is close to the WC 
state. Note here that such decrease was not observed in 
Ref. @, which investigated the quarter-filled system and 



reported quadratic increase of the relative fluctuations 
with V/t at strong interactions. In general stronger in- 
teractions are required for systems with smaller filling 
factors to form WC states, and accordingly it is expected 
that relative fluctuations begin to decrease at larger val- 
ues of V/t. We thus believe that further increase of the 
interaction strength should yield reduction of the relative 
fluctuations even in the quarter-filled case. 

Finally, we check the validity of the HF approxima- 
tion by comparing the HF results with those obtained 
via exact diagonalization of the Hamiltonian in Eq. (|^) 
for small systems. In Fig. [)](a) and (b), which displays 
the statistics of A 2 obtained via the two methods in the 
system of L = 4, one can see that the HF results are 
in good agreement with the exact ones except for some 
parameter regions: For V/t > 1, observed in Fig. ||(a) is 
that the HF approximation overestimates (A 2 ) for small 
W/V (< 0.3). Figure ^|(b) also shows such overestimation 
of the fluctuations in the HF approximation for the inter- 
mediate disorder strength 0.1 < W/V < 1 and V/t > 1. 
The distribution function plotted in Fig. g(c) and (d) pro- 
vides a clue to the overestimation: For W/V = 0.2 cor- 
responding to moderate disorder, the HF approximation 
generates more frequently large values of A 2 and hardly 
changes the minimum of A 2 , resulting in the overesti- 
mation of both (A 2 ) and <5A 2 . On the other hand, for 
smaller W/V = 0.1 the overall distribution merely shifts 
to the right with the shape almost unchanged. Accord- 
ingly, whereas the HF approximation tends to give larger 
values of (A 2 ), rather accurate values are obtained for 
8A2 in the case of small W/V, as shown in Fig. 
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FIG. 6: (a) Average, (b) fluctuations, and (c) relative fluc- 
tuations of the inverse compressibility as functions of W/V 
in the quantum system of size L — 8. The data marked by 
squares, circles, triangles, inverted triangles, and diamonds 
correspond to the hopping strength t/V = 0.1, 0.5, 1, 2, and 
10, respectively, while the thick lines represent the data in the 
classical limit (t/V = 0). 



IV. EFFECTS OF MAGNETIC FIELDS 

In this section we examine the effects of magnetic fields 
on the distribution and fluctuations of the inverse com- 
pressibility. Electrons hopping in applied magnetic fields 
acquire additional phases, leading to the kinetic term in 



the Hamiltonian 

H K = tY,( e ~ lA ' Jc h + e-^cXa). (7) 

The phase Ay associated with the hopping between sites 
i and j is given by 



A, 



±2irfyi for j = i ± x 
for j = i ± y 



in the Landau gauge, where / = Ba 2 /$>o is the frustra- 
tion parameter in the presence of a uniform transverse 
magnetic field B = Bi. We also adopt the free boundary 
conditions, which are convenient for describing the sys- 
tem with arbitrary values of the frustration parameter 
(i.e., under arbitrary magnetic fields) without any mis- 
match on the boundaries. Periodic boundary conditions 
in general cannot avoid such mismatch except for rather 
low-lying rational values of the frustration restricted by 
the system size. 

Some experimental groupstra have measured the con- 
ductance peak spacing in the presence of finite mag- 
netic fields, reporting that the distribution of A2 re- 
mains symmetric and Gaussian-like, as in the case 
of zero magnetic field, with slightly smaller width: 
(SA 2 / (A 2 ))b=o(«SA 2 / (A 2 )) b ^ - 1.2 ± 0.lj3 Note that 
those experiments probed the weak-field regime: The 
magnetic flux per electron in units of the flux quantum 
roughly amounts to 0.08^ or 0.002 to 0.04.0 

Motivated by this, we first concentrate on the weakly 
frustrated systems, where magnetic fields are weak 
enough not to change significantly the characteristics 
around the Fermi levels of corresponding noninteracting 
clean systems. Such a case is given by the frustration 
range < / < fdi, where fdi is the smallest nonzero 
value of the degenerate frustration defined below. We 
then consider the systems with large values of the frustra- 
tion, especially the fully frustrated system with / = 1/2. 

For clear identification of the magnetic-field effects, 
we first investigate the system of free electrons on a 
lattice, in the presence of a magnetic field. It is well 
known that in two dimensions the system is described 
by Harper's equation, displaying^the characteristic spec- 
trum of a very complex pattern. E3 In Fig. ^ we show the 
single-particle energy spectrum and the corresponding in- 
verse compressibility in the half-filling case. It is observed 
that the Fermi level becomes degenerate at some values 
of the frustration, leading the inverse compressibility to 
vanish. As the system size grows, the number of such 
degenerate frustrations, denoted by fd, increases rapidly 
and they spread out over the whole interval [0, 1/2) ex- 
cept near the full frustration (/ = 1/2). 

A recent studyi has revealed that in the absence of 
the magnetic field the degeneracy at the Fermi level in 
a free-electron system affects the inverse compressibility 
and its fluctuations. The single-particle energies in the 
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FIG. 7: Distributions of the normalized inverse compressibility in the quantum system of L — 8 for W/t = 1 and V/t = (a) 
0.1, (b) 1, (c) 10; for W/t = 5 and V/t = (d) 0.1, (e) 1, (f) 10. The solid and the dashed curves represent the best fit to the 
normalized Gaussian and to the WD distributions, respectively. 
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FIG. 8: Relative fluctuations of the inverse compressibility 
versus V/t in the quantum system of L = 8 for various values 
of W/t = 0.1(D), 0.5(B), 0.75(c), 1(.), 1.5(A), 2(a), 4(v), 
6(T), and 10(o). 



system are given by 



2t cos • 



L + l 



L + l 



(8) 



with n and m being integers from 1 to L. Accordingly, 
there are L degenerate states with (n, m) = (1, L), (2, L — 
1), . . . , (L, 1) at the Fermi level. This degeneracy yields 
finite average inverse compressibility together with non- 
vanishing relative fluctuations in the-.free-electron limit 
of the interacting disordered systemJ3 Such observation 
leads to a natural question whether the same phenom- 
ena - relatively small inverse compressibility and non- 



vanishing relative fluctuations - remain at nonzero de- 
generate values of the frustration and how the distribu- 
tion of the inverse compressibility changes with the mag- 
netic field. An answer to such questions are given below. 
In general, for large systems where the set of the de- 
generate frustration becomes very dense, the statistics of 
the inverse compressibility is expected to remain qualita- 
tively unchanged irrespectively of the frustration, except 
around / = 1/2. 



A. Weakly Frustrated System 

We now examine the effects of weak magnetic fields 
on the statistics of the inverse compressibility. Recalling 
that the influence of the magnetic field enters only in the 
hopping term as an additional phase, we can infer that 
the magnetic-field effects are not prominent in the sys- 
tem with small hopping strength. Indeed, for W/t > 2 
or V/t > 2, no distinctive change has been found numer- 
ically, at least in the presence of weak magnetic fields. 

Figure [ll] shows the dependence of (A2) /t on the frus- 
tration: The magnetic field initially enhances (A2) /t, 
which is expected from the increase of A2 with / in 
the free-electron system. As / approaches fdi, the first 
nonzero degenerate frustration in the corresponding non- 
interacting clean system, however, (A2) /t decreases to a 
value which is close to the zero-frustration value. Level 
repulsion due to random disorder smoothes out the de- 
pendence on /, especially at the degenerate frustration. 

On the other hand, fluctuations of the inverse com- 
pressibility show different behaviors with the frustration, 
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FIG. 9: Comparison of the HF results (filled symbols) and the exact ones (empty symbols with lines) for (a) the average and 
(b) fluctuations of the inverse compressibility in the quantum system of L = 4 for various values of V/t = 0.5, 1, 1.5, 2, 5, 
and 10 from below; comparison of the inverse-compressibility distribution for V/t = 10 and W/V = (c) 0.2 and (d) 0.1. The 
discrepancy between the two results is considerable in (d). 



depending on the strength of the interaction and disor- 
der. In the system with weak disorder [see Fig. |l2|(a)- 
(c) for W/t = 0.5], SA 2 /t versus / exhibits two minima 
for weak interactions V/t(< 2), which turn into a single 
minimum as the interaction is increased; these features, 
however, become less pronounced with the increase of 
the syst em size. For strong disorder (W/t — 2) shown in 
Fig. [13(d), (e), and (f), the peculiar structure observed 
in the weak-disorder regime fades away and 5A%/t tends 
to decrease with the frustration, particularly in large sys- 
tems. 

From the distribution of A2 presented in Fig. ^3|, one 
can get a hint as to the origin of the behavior of fluc- 
tuations in the weakly disordered system. Turning on 
the magnetic fields initially suppresses the appearance 
of small A2 and thus reduces fluctuations regardless of 
the interaction strength, as can be observed in Fig. |l3|(a) 
and (c). Further increase of /, however, raises signifi- 
cantly the portion of large A2 in the presence of weak 
interactions, inducing large fluctuations [see Fig. |i~3|(b)]. 
When / exceeds the value at which (A2) ft reaches its 
maximum, the reverse process takes place, producing an- 
other minimum of hA 2 /t. In contrast, strong interactions 
seem to suppress such extension of the distribution to the 



large-A2 region, leading to a curve with one minimum for 
SA 2 /t. 

To examine the relation between different responses to 
the magnetic field and localization properties, we intro- 
duce the participation ratio /3 and the inverse participa- 
tion number (IPN) P _1 (j/): 

p-» = j>Mr, (10) 

i 

where N e is the total number of electrons, pi is the expec- 
tation value of the electron density at site i, and ^i(v) i s 
the probability amplitude for the single-particle state \v) 
at site i. The participation ratio (3 is equal to unity in the 
state with a uniform density profile and decreases as the 
electron configuration becomes modulated or localized, 
eventually reaching the minimal value of the filling fac- 
tor N e /L 2 in the. limit that all the electrons are located 
at the N e sitesO The IPN, on the other hand, assumes 
unity for the level localized at one site and decreases as 
the state becomes delocalized, approaching its minimal 
value L~ 2 in the opposite limit of the uniformly extended 
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FIG. 10: (a) Single-particle energy spectrum and (b) the in- 
verse compressibility versus frustration in the half-filled non- 
interacting clean system of size L = 8. Since the spectrum 
possesses reflection symmetry with respect to the E = line, 
only lower half of the spectrum is shown here. 



stated 

We use these measures to characterize localization or 
modulation of the half-filled state \M) and search for the 
possible relation with the inverse compressibility, which 
reveals that for given frustration the participation ratio 
in general has negative correlations with the inverse com- 
pressibility: The smaller the inverse compressibility, the 
larger the participation ratio and vice versa. This obser- 
vation indicates that the electron profiles in the samples 
with smaller values of A2 tend to be uniform and thus 
explains the increase of A2 in the samples displaying rel- 
atively small A2 since such uniform profiles should be 
sensitive to the change in frustration. Localized states, 
on the other hand, are rather insensitive to frustration, 
and accordingly, the inverse compressibility in a sample 
with large A2 remains unchanged. To be more definite 
in the electron state, we have also calculated the IPN 
at the Fermi level and observed positive correlations be- 
tween the IPN and the inverse compressibility, thus con- 
firming the above argument. In the case that hopping is 
dominant over the interaction and disorder, the samples 
mostly have large values of the participation ratio and 
low values of the IPN at the Fermi level, implying sensi- 
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FIG. 11: Average inverse compressibility versus frustration / 
in the range [0, fdi] for various values of V/t =0(D), 0.5(o), 
1(A), 1.5(V), and 2(o) in the system of size L = (a) 4, (b) 6, 
and (c) 8. The disorder strength W/t is set to be 0.5. 



tivity to frustration. Variations of the sensitivity over the 
samples lead to the increase of width of the distribution. 

As / approaches fdi, the distribution reverses its be- 
havior, restoring that at / » 0. However, it should be 
noted that the distributions at / = and at / = fdi 
are not exactly the same. Rather, the distribution at 
/ = fdi is similar to the one at / = 0.005. Such dif- 
ference between the distributions at different degenerate 
frustrations may be attributed to the difference of the 
degree of degeneracy: At / = the system has the L- 
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FIG. 12: Fluctuations of the inverse compressibility versus frustration for various values of V/t =0(D), 0.5(o), 1(A), 1.5(V), 
and 2(o) in the weakly disordered system (W/t = 0.5) of size L = (a) 4, (b) 6, (c) 8 and in the strongly disordered system 
(W/t = 2) of size L = (d) 4, (e) 6, (f) 8. 



fold degeneracy at half-nllingQ but at / = fdi two-fold 
degeneracy is present, as can be seen in Fig. [l0](a). Since 
large degeneracy favors a smaller value of A2 even in the 
presence of the interaction and disorder, the probability 
for small A2 at / = is higher than that at / = fdi, 
which is verified in numerical calculations. 

To get more insight for large fluctuations at the degen- 
erate frustration in relation to the localization property, 
we have computed the average participation ratio (/?), 
which is plotted as a function of / for various interac- 
tion strengths in Fig. [14. The conspicuous decrease in 



((3) near the degenerate frustration signals an increase in 
spatial modulation or localization of the total electron 
density. Indeed in a perfectly clean system the partici- 
pation ratio is exactly equal to unity, irrespective of the 
frustration. As disorder is introduced, {(3) diminishes 
faster at the degenerate frustration. Such spatial local- 
ization at the degenerate frustration or derealization at 
non-degenerate frustration are believed to be related with 
the pattern of interference of electrons scattered at im- 
purities: Moving around the lattice under the magnetic 
field, the electrons acquire additional phases. At the de- 
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FIG. 13: Distribution of the inverse compressibility in the 
system of size L = 8 for W/t = 0.5. The interaction strength 
is given by V/t — 0.5 in (a) and (b) and V/t = 2 in (c); the 
frustration / is varied from to fdi- 



FIG. 14: Average participation ratio versus frustration in the 
range [0, /<n]. The system has the size L = (a) 4, (b) 6, 
(c) 8 together with the interaction V/t = 0(D), 0.5(o), 1(A), 
1.5(V), and 2(o). The disorder strength is given by W/t = 0.5. 



generate frustration such phases are expected to enhance 
the constructive interference. In connection with the de- 
generacy we can imagine linear combinations of degen- 
erate states at the Fermi level; any combination will be 
an eigenstate of the clean system, and among them we 
can choose more modulated or localized one which leads 
to less energy in the presence of disorder. At the non- 
degenerate frustration such choice is not possible, yield- 
ing rather extended states. 

The sensitivity of electron states to disorder at the de- 



generate frustration in turn generates larger fluctuations 
of the inverse compressibility over the samples. In Ref. ^ 
the enhanced fluctuations of the inverse compressibility 
were also attributed to localization of the wavefunctions 
around the edge of the dot. Here we stress that our argu- 
ment is not restricted merely to the half-filling case which 
gives fairly high degeneracy at / = 0. Figure 10(a) shows 
that there also exist many degenerate states at / = at 
other fillings. 

We conclude this section with some remarks on the 
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FIG. 15: (a) Average and (b) fluctuations of the inverse com- 
pressibility versus W/t at f — 1/2 (empty symbols with lines) 
and at / = (filled symbols) in the system of size L = 8 for 
various values of V/t. The data marked by squares, circles, 
triangles, inverted triangles, diamonds, and pentagons corre- 
spond to V/t — 0, 0.5, 1, 2, 5, and 10, respectively. 



comparison with the existing experiment. For the inter- 
action and disorder which are sufficiently strong but not 
so strong to overwhelm quantum effects, the fluctuations 
are observed to diminish with the frustration /, provided 
that / <C fdi- This leads to similar decrease in relative 
fluctuations, even for weak disorder and interaction, ow- 
ing to the increase in (A2). This result agrees quite well 
with the experimental results mentioned in the beginning 
of this section. Furthermore, the structure of the inverse- 
compressibility distribution has been found not to change 
appreciably, which also coincides with the experimental 
observation. 



B. Fully Frustrated System 

In the previous section, we have examined the effects 
of weak magnetic fields, up to the first degenerate frus- 
tration fdi, on fluctuations of the inverse compressibility. 
The arguments in the previous section as to the fluctua- 
tions at degenerate and non-degenerate frustrations seem 
to be applicable beyond fdi as well. However, as the sys- 



tem size grows, the degenerate frustrations proliferate in 
number to spread over the entire region of the frustra- 
tion and the energy splitting by the magnetic field at the 
Fermi level also shrinks. Accordingly, no distinctive dif- 
ference is expected to appear in the distribution of A2 at 
degenerate and non-degenerate frustrations. A remark- 
able exception happens at the full frustration (/ = 1/2), 
where quite large energy spacing is induced at the Fermi 
level in the noninteracting clean system (see Fig. |l^). For 
even L, the energy spacing reads 



iV2t 



cos 



ttL 



2(£ + l) 



(11) 



which decreases to zero in the thermodynamic limit 
(L — ► 00). Nevertheless, for moderate L the finite spac- 
ing given by Eq. ( |TT| ) discerns the full frustration from 
other frustration values, yielding qualitative difference in 
the distribution even for strong interactions. 

The average inverse compressibility and its fluctuations 
are displayed in Fig. [l5|. For V/t — 2 and 5, one can no- 
tice the disappearance of both the lumps in (A2) jt and 
the overshoots in 8/^.%/t, which develop at weak disor- 
der and strong interactions in the unfrustrated system 
(/ = 0). In Sees. II and [II, both the lumps in (A2) 
and the overshoots in 5A2 have been found to have their 
origin in the WC states due to the Coulomb interaction. 
Accordingly, the disappearance of such lumps and over- 
shoots implies that strong magnetic fields suppress the in- 
teraction effects on the electron distribution. To look into 
the distribution of the electrons, we have compared the 
average participation ratio at / = and that at / = 1/2 
in Fig. 16. It is observed that the latter, (/3)^ =1 / 2 , is 
always larger than the former, (/3)f_ , with the differ- 
ence significant for 2 < V/t < 5. Furthermore, (/3)y =1 y 2 
does not vary much with the interaction, up to V/t « 2; 
even for V/t = 5, its values are not so different from the 
values for V/t — 0. We thus conclude that the full frus- 
tration restrains the modulation of electrons induced by 
the repulsive Coulomb interaction between electrons and 
suppresses its influence on the fluctuations of the inverse 
compressibility. 

To verify the suppression of the Coulomb interaction 
beyond the HF approximation, we compute the partic- 
ipation ratio as a function of the frustration parameter 
via the exact diagonalization method in the clean system 
with a single impurity of strength W/t — 1, located at a 
corner of the 4x4 lattice, and plot the results for various 
values of V/t in Fig. 17, Note that small values of the 
participation ratio indicate the formation of the WC-like 
states due to interactions. We have intentionally inserted 
an impurity in the system in order to eliminate the de- 
generacy, which allows the two (degenerate) WC states 
to superpose and may keep the participation ratio still 
larger even when the electrons form a Wigner crystal. 
At / = and other degenerate frustration values, the 
average participation ratio is clearly shown to decrease 
with the interaction; even rather weak interactions lead 
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FIG. 16: Average participation ratio versus W/t at / = 1/2 
(empty symbols with lines) and f — (filled symbols) in the 
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FIG. 17: Average participation ratio versus the frustration / 
in the system of size L = 4 with a single impurity of strength 
W/t = 1, located at a corner of the lattice, for various values 
of V/t = 0.1(D), 0.5(H), l(o), 2(«), 5(A), 10(a), 20(V), and 
50(T). 



ubiquity of large values of A2 at / = 1/2 comes not from 
the Coulomb charging energy but from the large single- 
particle energy splitting due to the full frustration. On 
the other hand, for V/t = 0.5 and W/t = 2, Fig. ||(b) 
shows a wider distribution at / = 1/2, which may be 
attributed to the fluctuations of the degenerate states 
just above and below the Fermi level: The energies of 
the degenerate states fluctuate severely according to the 
disorder configurations, resulting in large fluctuations of 
the inverse compressibility. 

We conclude this section with a comment on the va- 
lidity of the HF approximation at / = 1/2. As observed 
in Fig. [It], the effects of the full frustration survive even 
in the presence of rather strong interactions, with the 
strength up to V/t = 10. Such effects are, however, not 
disclosed in the average and fluctuations of the inverse 
compressibility calculated via the HF method (see Fig. [l5| 
for V/t — 10). Further, the average participation ratio 
obtained within the HF approximation in the same sys- 
tem as in Fig. [l7| turns out to display considerable dis- 
crepancy at V/t — 10, which strongly suggests that the 
HF approximation breaks down for strong interactions 
V/t > 10. 



SUMMARY 



to appreciable reduction in the participation ratio. At 
the full frustration, in contrast, the participation ratio 
does not reduce much even in the presence of substantial 
interactions, indicating that the fully frustrated system 
is reluctant to be in a WC state and thus confirming the 
above HF results. 

Figure |l^ compares the distribution of A 2 at / = 
and that at / = 1/2. For V/t = 2 and W/t = 0.5, shown 
in Fig. 18(a), the width of the distribution at / = 1/2 
is much reduced compared with that at / = 0. Here the 



We have studied numerically the inverse compress- 
ibility and its fluctuations in two-dimensional Coulomb 
glasses, with emphasis on the quantum effects associ- 
ated with electron hopping as well as the frustration ef- 
fects due to applied magnetic fields. For a systematic 
study of this problem, we have begun with the corre- 
sponding classical system, revealing the detailed depen- 
dence of the inverse-compressibility fluctuations on the 
interaction and random disorder. Then the effects of 
electron hopping on the distribution of the inverse com- 
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FIG. 18: Distribution of the inverse compressibility at / = 
and / = 1/2 for (a) V/t = 2 and W/t = 0.5; (b) V/t = 0.5 
and W/t = 2. 

pressibility have been examined, particularly with regard 
to the interplay with the interaction and disorder. The 
results, obtained mostly via the Hartree-Fock approxi- 



mation, have been compared with those of existing ex- 
periments as well as of previous numerical studies: In 
addition to the symmetric Gaussian distribution with 
non-Gaussian tails, we have observed, in a sufficiently 
clean sample with strong interactions, a peculiar right- 
biased distribution of the inverse compressibility, which 
was indeed reported in experiment under similar situa- 
tions. Realization of even cleaner samples, where such 
an unusual distribution may be attainable with rather 
weak interactions, is thus expected to confirm the valid- 
ity of our results and to clarify the role of interactions. 
Remarkably, the relative fluctuations have been found 
to decrease eventually as the interaction strength is in- 
creased. This is contrary to the quadratic interaction 
dependence of the fluctuations, suggested in a previous 
numerical study. 

We have next investigated the effects of magnetic fields 
on the distribution and fluctuations of the inverse com- 
pressibility. Weak frustration generally suppresses the 
fluctuations; in a clean sample with weak interactions, 
however, fluctuations are enhanced in some range of the 
frustration. The difference in responses to the magnetic 
fields has been discussed in relation to localization prop- 
erties, which provides a good explanation for the effects 
of the magnetic fields on the distribution of the inverse 
compressibility. Finally, strong frustration has been ob- 
served to suppress the effects of the interaction. 
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